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Abstract 

The flow fields of unheated, supersonic free jets 
from convergent and convergent-divergent 
nozzles operating at M =0.99, 1.4 and 1.6 were 
measured using spectrally resolved Rayleigh 
scattering technique. The axial component of 
velocity and temperature data as well as density 
data obtained from a previous experiment are 
presented in a systematic way with the goal of 
producing a data base useful for validating 
computational fluid dynamics codes. The 
Rayleigh scattering process from air molecules 
provides a fundamental means of measuring flow 
properties in a non-intrusive, particle free manner. 
In the spectrally resolved application, laser light 
scattered by the air molecules is collected and 
analyzed using a Fabry-Perot interferometer (FPI). 
The difference between the incident laser 
frequency and the peak of the Rayleigh spectrum 
provides a measure of gas velocity. The 
temperature is measured from the spectral 
broadening caused by the random thermal motion 
and density is measured from the total light 
intensity. The present point measurement 
technique uses a CW laser, a scanning Fabry-Perot 
interferometer (FPI) and photon counting 
electronics. The 1 mm long probe volume is 
moved from point to point to survey the flow 
fields. Additional arrangements were made to 
remove particles from the main as well as the 
entrained flow and to isolate FPI from the high 
sound and vibration levels produced by the 
supersonic jets. In general, velocity is measured 
within ±10 m/s accuracy and temperature within 
±10°K accuracy. 


I. INTRODUCTION 

With the advent of the High Speed Civil 
Transport program, the flow field of the 
supersonic free jets has generated renewed 
interest. The sound generated from such jets is of 
primary importance. In contrast to a vast amount 
of acoustic information gathered over the last 4 
decades, there exists very little information on the 
flow field of such jets. The present work is a part 
of a collaborative program to predict supersonic 
jet noise using advanced computational 
techniques. The goal of the experimental part is to 
provide a reliable and accurate database, which 
can be used to validate the computational fluid 
dynamics and aeroacoustics codes. 

The lack of experimental data is primarily 
due to the problems in using the traditional 
experimental tools which are either too intrusive 
or use seed particles. Moreover, unlike its 
incompressible counterpart, the compressibility 
effects present in a high Mach number jet require 
a description of the thermodynamic quantities, 
such as density and temperature in addition to 
velocity. Simultaneous measurement of fluid 
dynamic and thermodynamic quantities is all but 
impossible using the traditional probes. The 
intrusive probes, such as the static and total 
pressure tubes or a hot-wire probe cause 
significant alterations that can be easily detected 
as a change in the noise emitted by the jet. The 
optical seed particle based techniques, such as the 
Laser-Doppler velocimetry (LDV) and Particle 
Image Velocimetry (PIV) are partially successful 
in providing time averaged velocity data . 


NASA/TM— 2004-212391 


1 


However, various biasing errors encountered in 
turbulent flows, problems associated with seed 
particles passing through shock waves, beam 
steering effects from sharp density gradients and 
problems in measuring high Doppler frequency 
shift using commercially available electronics 
make such techniques either difficult or of 
doubtful value in supersonic flows. The Rayleigh 
scattering technique, on the other hand, depends 
on the light scattering from air molecules. Since, 
no seed particles are used and only a single beam 
is necessary, many problems associated with LDV 
and PIV can be overcome. Moreover, information 
is obtained directly from the molecular 
distributions; therefore, the Rayleigh scattering 
technique is a more fundamental way of 
measuring flow properties. 

In the last few years, a number of Rayleigh 
scattering based optical flow diagnostic 
techniques had been developed. To measure flow 
velocity and temperature one needs to spectrally 
resolve the Rayleigh scattered light, which has 
been pursued by two means. The first method is to 
use a tunable laser and molecular absorption 
filters, such as the hyperfine absorption lines of 
Iodine vapors (Forkey, Lempert & Miles 4 , Elliot 
& Samimy 5 , Seasholtz, Buggele & Reeder 6 and 
others). The second method is to use a Fabry- 
Perot interferometer (Seasholtz, Zupanc & 
Schneider , Snyder et al . Pitz, Cattolica & 
Talbot 9 , Lock, Seasholtz & John 10 and others). 

The two methods have their own advantages and 
disadvantages. The present work follows the 
second methodology, which perhaps is better 
suited to obtain unsteady time accurate 
measurements, a long-term goal of the current 
program. A combination of narrow line width 
continuous wave laser, scanning Fabry-Perot 
interferometer and photon counting electronics is 
used to perform point measurement of velocity 
and temperature in supersonic flows. It should be 
mentioned here that a simpler version of the setup 
(without the Fabry-Perot interferometer) has been 
successfully used to measure density in similar 
flow fields. This has been reported earlier (Panda 
& Seasholtz 11 ). 


Since it is relatively easy to produce a 
small supersonic jet in the laboratory 
environment, free jets have been common test 
grounds for various optical techniques including a 
few variations of the Rayleigh scattering based 
techniques 4 ' 6 . The velocity and density data 
reported in the above studies are from a few 
points, and are insufficient to characterize the 
global flow field, which is the motivation of the 
current study. The primary attention of the earlier 
work was a demonstration of the proof of 
principle where the optical setup is at primary 
focus and the purpose of a tiny free jet is to create 
a supersonic flow. In the present work additional 
complexities, such as high noise level, vibration, 
particle contamination, condensation etc. present 
in a comparatively larger facility are addressed. 

In addition to the velocity and temperature 
measurements density data obtained from the 
same flowfield are also presented in this paper for 
the sake of completeness. Although it is 
straightforward to measure all three quantities 
simultaneously, the current setup was used to 
measure velocity and temperature only and 
density values were measured separately 11 ' 12 . For 
density data, intensity of the total scattered light 
needs to be measured. Since the Fabry-Perot 
interferometer is unnecessary for this purpose, a 
higher accuracy level is achievable in separate 
measurement. 

II. THEORY OF THE 
MEASUREMENT TECHNIQUE 

Rayleigh scattering is a consequence of elastic 
light scattering process from gas molecules. 
Traditionally, molecular scattering has been 
separated into two regimes, kinetic and 
hydrodynamic (also known as Brillouin 
scattering) based on the relative importance of 
scattering from the very high frequency density 
fluctuations with a characteristic length scale 
comparable to the molecular free path. For the 
present experimental conditions, involving not 
very high temperature and pressure and 
approximately 90° scattering, the influence of the 
Brillouin scattering is relatively small and the 
scattering process can be approximated as a sum 
of contributions from individual molecules present 
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at the probe volume. In this limit, the light 
scattered from each of the molecules is Doppler 
shifted according to its velocity and direction. 

A simplistic description of the 
measurement process using laser induced 
Rayleigh scattering is schematically shown in Fig. 

1 . Since the Doppler shift frequency is of the order 
of a few GigaHertz, a narrow line width (of the 
order of 10MHz) incident laser beam is necessary 
to resolve Rayleigh spectrum. In the case of a 
moving gas media, the bulk motion is 
superimposed on the random velocity of the 
individual molecules; therefore separation 
between the peaks of the incident laser line and 
the Rayleigh spectrum provides a measure of the 
bulk velocity. The Full Width at Half Maxima 
(FWHM) depends on the spread of molecular 
motion and, therefore, is a measure of gas 
temperature. Finally, the total light under the 
Rayleigh spectrum is proportional to the 
molecular number density and provides a measure 
of gas density. Therefore, a single Rayleigh 
spectrum carries infonnation of one component of 
gas velocity, bulk temperature and bulk density. 

Figure 2 presents a scattering diagram, 
where kj is the incident wave vector, k s is the 
scattered wave vector pointing towards the 
collection optics, 9 S is the scattering angle and, k 
= k s - kj, is the scattering vector. For k s « kj, the 
magnitude of the scattering vector becomes k = 2 
kj sin(0 s /2), and the Doppler shift from the bulk 
velocity U of the flow is 

. TT 2 U cos® sin(— ) 
f d = — = 2_ (!) 

2 n Ai 

where, ® is the angle between the scattering 
vector k and the velocity vector U, and k; is the 
wavelength of the incident light. The same 
Doppler shift principle applies to the individual 
molecules as well. The gas molecules are assumed 
to have a Maxwellian velocity distribution f(V): 

f(V;d 3 V = ^Aexp(-^^-)d 3 V, (2) 

n a a 

where, n = the molecular number density, a = 
"most probable speed" = (2k B T/m) 1/2 , k B is the 


Boltzmann's constant, m is molecular mass and T 
is gas temperature. If Brillouin scattering has 
negligible contribution, then the bandwidth 
(FWHM) of Rayleigh spectrum, due to the 
Maxwellian molecular velocity distribution 
becomes: 

Af= 0.265 k(—^—) 1/2 . (3) 

m 

Clearly, for a fixed optical setup (k fixed) and a 
given gas mixture (m fixed) the bandwidth of 
Rayleigh spectrum is proportional to the square 
root of temperature. 

It should be pointed out here that the above 
description of velocity, temperature and density 
measurement is somewhat simplistic. The Tenti 
S6 model 13 used in data analysis considers both 
Rayleigh and Brillouin parts of scattering. Instead 
of just measuring the FWHM of Rayleigh spectra 
the Tenti's model function is fitted through 
experimental data to improve accuracy of 
temperature measurement. Since density depends 
on the total scattered light, no spectral analysis is 
required. 


Fabry-Perot Interferometer. The high resolution 
required in Rayleigh spectroscopy necessitates the 
use of a Fabry-Perot interferometer 1415 (FPI). It 
acts as a light filter whose transmittance is 
dependent on the frequency f of the incident light, 
the angle q of the incident light with the optical 
axis, the refractive index p in the medium between 
two flat mirrors and the mirror separation d (Fig. 
3). In an ideal no loss situation the transmittance is 
given by the Airy function: 


I. = [l + (- 


2N 


>w< 2*— r , ^ (4) 


n A { 

where, N is the interferometer finesse that 
detennines the sharpness of the filter roll-off; a 
higher finesse translates into a sharper roll-off. For 
an ideal interferometer finesse is a function of the 
mirror reflectance only, however; in reality factors 
such as imperfections in the mirror coating, 
misalignment of the optical train, vibration and 
noise (both acoustic and electronic) degrades the 
finesse. An incorrect finesse value affects the 
measurement accuracy and one needs to monitor 
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finesse continuously during the course of 
experiment. Figure 3 shows that the light coming 
off an optical fiber is imaged through the FPI 
using a focusing and a fringe forming lens. As 
light from different part of the fiber makes 
different 9 angle with the optical axis, even a 
single frequency light is imaged as concentric 
rings at the image plane. In the present setup, to 
measure frequency of incoming light, the FPI was 
operated in the scanning mode, that is, the mirror 
separation d is continuously changed by moving 
one of the two mirrors. A scan that changes the 
mirror spacing by k,/2 causes the FPI output to 
move from one transmission peak to the next; a 
separation represented by the free spectral range 
of the interferometer: 

c 

Free spectral range = — , (5) 

where, c is the speed of light. As the plates are 
scanned the fringes formed at the image plane 
move inward. In practice, additional spurious 
concentric rings are found outside the fiber face. 
Based upon the optimization calculation of 
Hernandez 15 the diameter of the pinhole, placed at 
the image plane, is chosen to be 30% of the first 
fringe diameter. 

The frequency spectrum S(v) measured 
using the Fabry-Perot interferometer is a 
convolution of the actual Rayleigh scattering 
spectrum Sr(u) and the characteristic of FPI filter 
represented as a Instrument function I(v): 

SmO) = Gr|I(v')S r (v- v')dv' +B, (6) 

where, Gr is the total Rayleigh scattered power 
and B is the uniform background from the 
photomultiplier dark current. The laser line width, 
finesse and free spectral range are included in the 
instrument function. Note that frequency is 
represented as v which is a relative scale referred 
to an arbitrary reference value. This is allowable 
since the frequency difference between the 
incident and the scattered light is of interest for 
aerodynamic measurements. The instrument 
function also includes the location of the 
transmission peak, relative to which the Doppler 
shift is measured. Since the instrument function 
tends to change with a small drift in the ambient 


condition and the instrument controlling 
electronics, it has to be monitored continuously. 
Note that, unlike density measurement, velocity 
and temperature measurement using Rayleigh 
scattering technique do not require any fluid 
dynamic calibration. However, the slow drift of 
FPI and the laser line requires reevaluation of the 
instrument function before each Rayleigh 
spectrum is acquired. 

III. EXPERIMENTAL SETUP 

The experiments were conducted at the NASA 
Glenn Research Center. Two convergent- 
divergent nozzles were designed using the method 
of characteristics for fully expanded operating 
Mach numbers of M = 1.4 and 1.8. Later on, 
centerline surveys showed that the least shock 
operation is achieved at respectively, M = 1.395 
and 1.795. The operating conditions are provided 
in Table I. The exit diameter for both nozzles is 
25.45mm (1.002 inches). Figure 4 shows a 
schematic of a CD nozzle block. The axial and 
radial coordinates of the inside contour of the C-D 
nozzles are provided in Table II. In addition a 
convergent nozzle of 25.4 mm exit diameter was 
used to produce M = 0.99 jet. Figure 5 shows a 
schematic of the Rayleigh scattering setup. An 
additional view showing a part of the arrangement 
around the nozzle block is shown in photograph 6. 
The supply air was unheated with an average total 
temperature of 24° Celsius. This was also the 
average ambient temperature. 

Optical configuration: In essence, the Rayleigh 
scattering setup consists of a continuous laser 
beam focused at a point (probe volume) in the 
flow field; the light scattered by the air molecules 
is collected, analyzed by a Fabry-Perot 
interferometer and measured by a photomultiplier 
tube. The extremely high noise level produced by 
the supersonic jets was a concern for the stable 
operation of the laser and the Fabry-Perot 
interferometer (FPI). The optical alignment of the 
flat mirrors present in the interferometer is 
especially sensitive to noise and vibration. Any 
vibration capable of producing a few nanometers 
of displacement in the mirror alignment creates a 
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measurable degradation in finesse. Since, a 
successful operation of the Rayleigh scattering 
system was critically dependent on the proper 
management of the noise and vibration issue, it 
was decided to split the optical arrangement into 3 
parts. The first part is the laser light source: a 
Neodymium Vanadate (Nd:YV04) laser that 
produced 5Watt power in 532 nanometer 
wavelength and has a narrow 100 MHz bandwidth 
was used. The laser source was kept in the same 
room that houses the jet facility. However to 
ensure stable operation an anechoic box was built 
around the laser and the light was passed through 
30 meter long, 0.365tnm core diameter, 
multimode, fiber-optic cable (Fig. 5). 

The second part consisted of an 
arrangement built around the jet facility and 
mounted on a two-axis traversing mechanism that 
facilitated point by point movement along the flow 
(x-axis) and transverse (r-axis) directions (Figs. 5, 
6). The transmitting side of the optical train was 
designed to focus the incident beam at the 
measurement volume using an f/2, 80tnm diameter 
and an f/3.3, 82tnm diameter achromatic lens. The 
scattered light was collected from a scattering 
angle 0 S = 85° (Fig. 2) using identical lenses and 
was coupled to a 0.55mm core diameter fiber. The 
optical configuration was designed to align the 
scattering vector along the jet flow direction so 
that the axial component of the velocity vector is 
measured directly. The f/3.3 collection lens fixed 
the collection aperture. The laser power delivered 
to the probe volume was 3.5 Watts. However, the 
laser light at the probe volume was randomly 
polarized by its transmission through the optical 
fiber. Thus, because of the polarization sensitivity 
of the Rayleigh scattering, only about 1/2 of the 
laser power was effectively used. The core 
diameter of the collection fiber and the 
magnification factor of the collecting lenses fixed 
the probe volume length to 1.03mm. The diameter 
of the beam waist was 0. 16mm. 

As mentioned earlier, the Fabry-Perot 
alignment and the laser operation require constant 
monitoring in the form of evaluating the 
instrument function I(v) before analyzing an 
Rayleigh spectrum. For this purpose a small 


amount of incident laser light coming off the 
transmitting fiber is coupled to the collection fiber 
using an additional reference fiber (Fig. 5). A 
small part of retro-reflection from the transmitting 
lenses is coupled at the input of the reference 
fiber. The output is mounted on a pneumatic 
actuator that can be moved to bring the fiber end 
in front of the receiving fiber. At this positron 
Rayleigh scattered light is blocked and the 
receiving fiber transmits a small amount of the 
incident laser light. This light is analyzed by the 
FPI to determine I(v). To couple Rayleigh 
scattered light the actuator is moved back. The 
actuator motion was slowed down using additional 
pneumatic control to minimize vibration and 
shock from quick deployment and retracing. The 
on-off motion of the actuator was controlled 
through a 5V TTL signal. To minimize stray light, 
the incident beam was directed into a beam dump. 
A similar arrangement was made in the collection 
side to provide a dark background. The jet 
operating conditions were monitored via a 
pressure transducer and a thermocouple probe 
placed in the jet plenum. 

The third part of the overall setup consisted of 
the interferometric arrangement for the spectral 
analysis of the collected light. This was kept in an 
isolated loft in an adjoining room that was free 
from the high noise level produced by the jet. The 
collected light was brought from the test cell using 
a fiber-optic cable. The Fabry-Perot interferometer 
used is a commercially available one with 70 mm 
diameter mirrors, k/200 flatness and 90% 
reflectivity. The mirrors are set at 20.3 mm 
spacing, producing a free spectral range (FSR) of 
7.4 Ghz. The frequency range scanned was about 
2 FSR which correspond to a change in the mirror 
spacing of X. The light coming out of the 
collection fiber is collimated using a f/2, 80 mm 
diameter achromatic lens into the FPI, and an f/5 
100mm diameter fringe fonning lens was used to 
focus FPI output on the pinhole. Light from the 
pinhole is once again collimated and measured 
using a thermo-electrically cooled photomultiplier 
tube. 

A few other hurdles had to be overcome 
before the Rayleigh scattering technique could be 
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successfully implemented. The most difficult one 
arose from the dust particles. Since particle (Mie) 
scattering is a few orders of magnitude stronger 
than the Rayleigh scattered light, presence of dust 
particles can completely swamp the Rayleigh 
signature. The dry air, supplied to the rig from a 
central high pressure facility, was passed through 
two sets of filters for dust removal. This made the 
primary jet air very clean; however particles from 
the entrained air made density measurement 
impossible in the jet shear layer. To mitigate this 
problem, an air filter system that provided dust 
free air in the jet neighborhood was installed. The 
clean air was routed via a flexible duct and 
exhausted vertically at a very low speed (~ 8 m/s) 
through a 30.5cm X 38cm hood over the nozzle 
block (figs. 5, 6). The clean air from the hood 
flows in a vertical direction, nonnal to the main 
jet, and in effect, creates a crossflow situation. 
However, any perturbation to the supersonic jet is 
expected to be minimal due to the very low 
crossflow velocity. 

The electrical circuit supplying actuation 
voltage (ramp signal) for the piezo-electric 
elements in the FPI was found to be susceptible to 
electrical noise in the house power supply. A line 
filter circuitry has to be used to eliminate such 
noise. 

Data acquisition and processing: The low level 
of the scattered light called for a photon counting 
system and the cooled photomultiplier tube (PMT) 
was especially suitable. The PMT output was 
tenninated on a 50 Ohm terminator, amplified by 
a wide bandwidth electronic amplifier and the 
final pulse counting was perfonned through a 
photon counter. Typically, the FPI was set to scan 
two free spectral ranges over half a second time, 
which was divided into 100 smaller intervals 
(bins) and the number of photoelectrons arriving 
during each of the smaller intervals were counted. 
To reduce the effect of random noise in scanning 
and to improve photon count averaging over 4 
successive scans were perfonned. Data acquisition 
at each measurement station was a two step 
process. 


First a reference spectrum was obtained to 
detennine the instrument function. For this 
purpose the actuator was moved to position the 
output end of the reference fiber in front of the 
receiving fiber. A small amount of incident laser 
light was transmitted to the Fabry-Perot 
interferometer and a reference spectrum was 
collected (Fig. 7a). In the second step the 
pneumatic actuator was retracted back, allowing 
Rayleigh scattered light to be transmitted through 
the receiving fiber and a Rayleigh spectrum was 
collected (Fig. 7b). Such a two-step acquisition 
process was repeated for all data points while the 
probe volume was moved from point to point in 
the flow field. 

The post-processing of the acquired data is 
also a two step process. First a model Airy 
function (similar to that of equation 4) is fitted 
through the reference signal using a maximum 
likelihood curve fitting procedure to obtain the 
instrument function. Next the Rayleigh spectrum 
is fitted to a model function represented by 
equation (6) using the same curve fitting 
procedure where Poisson's statistics are assumed 
for the photon arrival rate. For the Rayleigh 
spectrum Tenti's S6 model was used . The curve 
fitting procedure provides estimates of the 
parameters of interest: velocity, gas temperature, a 
relative measure of gas number density (Gr in 
equation 6), and uniform background (B in 
equation 6). 

IV. RESULTS AND DISCUSSION 

Validation: Figure 8 shows some measured data 
obtained from the core region (r/D = 0, x/D = 2.5, 
where r represents radial distance from the jet 
core, x is the downstream distance from the nozzle 
exit and D is the nozzle exit diameter) of subsonic 
free jets. The jet core is particularly suitable for 
validation as the flow properties remain uniform. 
The jet Mach number M was changed by changing 
the plenum pressure Po. The axial velocity u, 
temperature t and density r at the core of a free jet 
emitting from an convergent nozzle can be 
detennined easily using isentropic relations until a 
choked condition (M < 1) is reached: 
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where. To is the plenum temperature, P am b is the 
ambient pressure, R is the universal gas constant 
and g is the ratio of specific heats. 

The velocity and temperature data follows 
the expected values closely. At a fixed jet 
operating condition multiple Rayleigh spectra 
were obtained to determine repeatability which 
also provides an uncertainty estimate. It is found 
that the velocity data are in general repeatable 
within ±10 m/s. The measured temperature is 
found to be somewhat lower than the actual 
temperature; the bias is about -8°K, and the 
superimposed random uncertainty is additional 
±8°K. A percentage uncertainty estimate is 
perhaps inappropriate; as there is no systematic 
uncertainty, that is, the random error is found to 
be independent of the measured velocity or 
temperature. The primary sources of uncertainty 
are traced to the variations in incident laser 
frequency and FPI mirror scanning using piezo- 
electric transducers (PZT). The shot noise 
contribution is expected to be small due to the 
relatively high rate of photon arrival. Some of the 
data points presented later in this paper are 
affected by the presence of small amount of dust 
particles in the flow. The dust particles deteriorate 
temperature measurements while velocity data 
remain unaffected. Noise associated with the 
mirror scanning was found to degrade the 
instrument finesse and had contributed towards 
the uncertainty in fitting model functions. Finally, 
the uncertainty level associated with the density 
measurements is the lowest, ±1%, among all. 

Fig. 9 shows centerline decay of velocity, 
temperature and density measured in the M = 0.99 


jet. The density measurements are from a different 
Rayleigh scattering arrangement 1 1 ’ 1 2 . In the 
present unheated jets isentropic expansion through 
the nozzle cools the core flow below ambient and 
increases density above ambient condition. 

Mixing with the ambient air, however, causes a 
rise in temperature and fall in density, both of 
which approach the ambient condition. These 
features are seen in the presented data. The 
uncertainty in velocity and temperature 
measurement has given rise to the random data 
scatter. There are no shock cells in this jet as 
apparent from the density data that has the lowest 
uncertainty level. 

Fig. 10 presents radial surveys obtained 
from various axial stations in the same M = 0.99 
jet emitting from a convergent nozzle. All 3 
profiles - velocity, temperature and density - show 
top hat profiles at the first x-station at x/D=2 and 
progressive change towards the Gaussian profile. 
Among the 3 data sets, the temperature data are 
the noisiest. The uncertainty level increases 
significantly beyond about 7.5D downstream 
distance. The reason lies with an increase in the 
number of particles passing through the probe 
volume. It has been pointed out earlier that the 
clean air supply hood extends about 7.5D from the 
nozzle exit, after which particle laden room air 
progressively mixes with the main airflow. 
However, the particles are still few as reflected in 
the total scattered light used for density 
measurement. Kourous and Seasholtz 16 have 
demonstrated that velocity and temperature 
infonnation can be reliably extracted when the 
flow is laden with a large number of particles. 
When there are only occasional particles and a 
scanning Fabry-Perot is used to resolve the 
Rayleigh spectrum , such as in the present 
situation, uncertainty level in temperature 
measurement significantly increases. The reason 
lies with the ability to model the particle light. A 
large number of particles produce a continuous 
distribution that is spread around the peak in the 
Rayleigh spectrum and is separable from the 
molecular scattered light. In contrast, an 
occasional particle affects only the spectral part 
acquired at the instance of particle passing, and 
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therefore, a suitable modeling of the particle 
scattered light becomes impossible. The outcome 
is the increased level of uncertainty in 
temperature measurement as seen in the profile 
obtained at x/D = 10 in Fig. 10 (b). 

The centerline and radial surveys in the 
supersonic plumes of M =1.4 and 1.8 C-D nozzles 
are shown in the figures 1 1 through 14. Note that 
the data were obtained at the least shock operating 
conditions of respectively, M = 1.395 and 1.795. 
The least shock operating conditions were 
determined from the centerline density surveys 
(similar to Figs. 11 & 13) made at various 
pressure ratio conditions. Since, density 
measurement is most accurate among the three, 
periodic modulations from the shock cells can be 
clearly identified from this data set. It should be 
pointed put that the data presented in Figs. 1 1 
through 14 are from an earlier setup that used an 
Argon-ion laser. Occasional mode hops 
encountered in that laser operation has caused 
higher uncertainty level in velocity measurements 
in a few stray data points. Figs. 1 1 (c ) and 13(c ) 
show the presence of weak cell patterns, 
consisting of compression and expansion waves, 
in the initial region of the jet plumes. Such 
patterns are not clearly discernable in the velocity 
and temperature data due to the random 
uncertainty. The shock cell patterns are somewhat 
stronger in the M = 1.4 jet, where the radial 
profiles (Fig. 12) also show distortions due to the 
shock cells. In a perfectly expanded jet, the radial 
profiles close to the nozzle exit are expected to be 
top hat in shape; presence of any shock cell causes 
distortion of the uniform part. Noticeably random 
uncertainty in the temperature data, at the farthest 
downstream stations for both supersonic jet cases, 
are not as high as for the M=.99 jet. The reason 
lies with the slower spreading rate expected with 
an increase in the Mach number, which in turn 
reduces the number of entrained particles. 

V. SUMMARY 

A Rayleigh scattering setup using the green line 
from a CW Argon-ion laser, a scanning Fabry- 
Perot interferometer and photon counting 


electronics were setup to perform point 
measurements in high-speed jet flow fields. The 
fundamental principle and an outline of the 
measurement process are described. To alleviate 
problems associated with high noise levels emitted 
by the jets the laser light source was placed in an a 
anechoic box and the Fabry-Perot interferometer 
with the associated optics was moved to an 
adjoining test cell that provided a quiet 
environment. Light transmission was performed 
through fiber-optic cables. To avoid problems 
associated with particle scattering filtered air was 
used for both the primary flow and in the 
immediate neighborhood of the jet. From 
measurements performed under known flow 
conditions, uncertainties in velocity measurement 
is established as ±10 m/s. The measured 
temperature is found to be lower by 8°K with a 
superimposed random uncertainty of ±8° K. 

Two axisymmetric convergent-divergent 
(C-D) nozzles with operating Mach numbers 1 .4 
and 1.8 were designed and fabricated. The plumes 
from these two C-D nozzles and a third case 
involving a convergent nozzle operated at M=0.99 
were surveyed using the Rayleigh scattering 
technique. In addition to the velocity and 
temperature data earlier density measurements 
perfonned in the same plumes are presented in 
this paper. The centerline decay characteristics 
and the radial profiles obtained up to a 
downstream distance of 1 5 nozzle diameters are 
expected to be valuable for computational fluid 
dynamics code validation. 

References 

'figgins, P. L. & Jackson, D. A. 1974 Laser-Doppler 
velocity measurements in under-expanded free jet. 
J. Phys. D: Appl. Phys. 7, 1894-1906. 

‘Ross, C., Lourenco. L. M., and Krothapalli, A. 
1994 Particle image velocimetry measurements in a 
shock-containing supersonic flow. AIAA 94-0047. 

Wemet, M. P. 1995 Fuzzy inference enhanced 
infonnation recovery from digital PIV using cross- 
correlation combined with particle tracking. NASA 
TM 106896. 


NASA/TM— 2004-212391 



4 Forkey, J. N., Lempert, W. R. and Miles, R. B., 
"Accuracy limits for planer measurements of flow 
field velocity, temperature and pressure using 
filtered Rayleigh Scattering," Experiments in Fluids, 
Vol. 24, 1998, pp. 15 1-162. 

5 Elliott, G. S. & Samimy, M. 1996 Rayleigh 
scattering technique for simultaneous measurements 
of velocity and thermodynamic properties. AIAA J., 
34 (11), 2346-2352. 

6 Seasholtz, R. G., Buggele, A. E. and Reeder, M. 
F., "Flow measurements based on Rayleigh 
scattering and Fabry-Perot interferometer," Optics 
and Lasers in Engineering, Vol. 27, pp. 543-570. 
7 Seasholtz, R. G., Zupanc, F. J. & Schneider, S. J. 
1992 Spectrally resolved Rayleigh scattering 
diagnostic for hydrogen-oxygen rocket plume 
studies. J. Propulsion & Power, 8(5), 935-942. 

o 

Snyder, S. C., Reynolds, L. D., Lassahn, G. D., 
Fincke, J. R., Shaw, G. B., Jr. and Kearney, R. J., 
"Determination of gas-temperature and velocity 
profiles in an argon thermal-plasma jet by laser light 
scattering," Physical Review E, Vol. 47, no. 3, 1993, 
pp. 1997-2005’. 

9 Pitz, R. W., Cattolica, R. Robben, F. and Talbot, 
L., "Temperature and density in a Hydrogen- Air 
flame from Rayleigh scattering," Combustion and 
Flame, Vol. 27, no. 3, 1976, pp. 313-320. 

10 Lock, J. A., Seasholtz, R. G. & John, W. T. 1992 
Rayleigh-Brillouin scattering to detennine one- 
dimensional temperature and number density 
profiles of a gas flow jet. Appl. Opt., 31(15), 2839- 
2848. 

n Panda, J. and Seasholtz, R. G., " Density 

measurements in underexpanded supersonic jets 

using Rayleigh scattering," AIAA paper no. 98- 

0281,1998. 

12 

Panda, J. and Seasholtz, R. G., “ Density 
fluctuation measurement in supersonic fully 
expanded jets using Rayleigh scattering. AIAA 99- 
1870, 1999. 

l3 Tenti, G., Boley, C. D. and Desai, R. C., "On the 
kinetic model description of Rayleigh-Brillouin 
scattering from molecular gases," Canadian J. of 
Physics, vol. 52, no. 4, Feb. 1974, pp. 285-290. 


14 Vaughan, J. M., The Fabry-Perot Interferometer, 
History, Theory, Practice and Applications. Adam 
Hilger, Bristol, 1989. 

15 Hemandez, G., Fabry-Perot Interferometers, 
Cambridge University Press, Cambridge, UK 1986. 
16 Kourous, H. E. and Seasholtz, R. G., “Fabry-Perot 
interferometer measurement of static temperature 
and velocity for ASTOVL model tests, NASA TM 
107014, 1994. 


Table I. Operating conditions 


Sp. heat ratio, y 


1.4 


Total temp., T 0 


300°K 

Ambient temp., T a 


297°K 

Ambient density, p a 


1.16 Kg/m 3 

Ambient Pressure, P a 


98.95 KPa 


M=0.99 

Conv. 

Nozzle 

M=1.4 

nozzle 

M=1.8 

nozzle 

Minimum shock at M = 

- 

1.395 

1.795 

Plenum pressure, Po 

(KPa) 

185.1 

312.7 

564.2 

Jet density, p, (Kg/m 3 ) 

1.37 

1.6 

1.89 

Static temp, T, (°K) 

251 

216 

182.4 

Jet velocity, Uj (m/s) 

314 

411 

486 
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Laser 

line 



Relative frequency in GHz 

Fig. 1 . Schematic to explain principle of 
Rayleigh scattering technique. 
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TableP. Nozzle coordinates. 


M - 1.8 CD nozzle 

M * 1.4 CD nozzle 

x (inches) 

radius (r.inches) 

x (inches) 

r (inches) 

-5.0 

1.5 

-4.926 

„ 

-3.773 

1.338 

-4.366 

1.482 

-2.593 

0.962 

-3.905 

1.437 

-2.157 

0.752 

-3.444 

1.368 

-2.049 

0.691 

-2.938 

1.257 

-1.943 

0.624 

-2.487 

1.132 

-1.761 

0.511 

-2.082 

0.992 

-1.662 

0.47 

-1.602 

0.789 

-1.535 

0.436 

-1.356 

0.664 

-1.464 

0.425 

•1.148 

0.552 

-1.391 

0.4188 

-0.98 

0.499 

-1.361 

0.4181 

-0.82 

0.4763 

-1.327 

0.4183 

-0.786 

0.4751 

-1.297 

0.419 

-0.749 

0.4749 

-1.244 

0.422 

-0.735 

0.475 

-1.189 

0.428 

-0.72 

0.4755 

-1.161 

0.431 

-0.685 

0.4772 

-1.118 

0.436 

-0.66 

0.4786 

-1.076 

0.441 

-0.635 

0.4798 

-1.033 

0.446 

-0.575 

0.4821 

-0.883 

0.46 

-0.546 

0.4834 

-0.777 

0.47 

-0.489 

0.4863 

-0.663 

0.479 

-0.468 

0.4871 

-0.577 

0.485 

-0.379 

0.4895 

-0.505 

0.489 

-0.255 

0.4938 

-0.289 

0.498 

-0.145 

0.4982 

0.0 

0.501 

0.0 

0.501 



Fig. 4. Schematic of the nozzle block. Coordinates of the inside contour are provided in Table 


II. 
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Nd: Vanadate laser 
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Fig. 5. A schematic of the experimental arrangement. 
























Fig. 6. Photograph of arrangement around jet facility. 



Relative frequency v in GHz 

Fig. 7. Sample of (a) reference (b) Rayleigh scattering 
spectra (centerline M=0.99 jet); solid line: fitted model 
function, symbols: experimental data. 
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Fig. 8. A comparison between measured (symbols) 
and expected (solid line) variation in M = 0 to 1 range 
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Density (Kg/m 3) Temperature (K) Velocity (m/s) 







velocity (m/s) 



temperature (deg K) 



density (Kg/m* *3) 


Fig. 10. Radial profiles from indicated axial stations in M = 0.99 jet. The chained lines indicate 
ambient condition. Each profile is separated by the distance between 2 chained lines. 
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Density (Kg/m**3) Temperature (K) Velocity (m/s) 




Fig. 11. Centerline variations in M = 1 .395 jet from a C-D nozzle. 
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2 4 

density (Kg/m* *3) 


Fig. 12. Radial profiles from indicated axial stations in M = 1.395 jet. The chained lines indicate 
ambient condition. 
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Density (Kg/m**3) Temperature (K) Velocity (m/s) 



0 5 10 15 20 

x/D 

Fig. 13. Centerline variations in M = 1.8 jet 


Note (added after publication): The centerline temperature decay of Fig. 13(b) has significantly higher 
uncertainty due to slight condensation in the jet. All other data have uncertainty mentioned in the text. 
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0 500 


velocity (m/s) 



2 

density (Kg/m**3) 


Fig. 14. Radial profiles from indicated axial stations in M = 1.8 jet. The chained lines indicate 
ambient condition. 

Note (added after publication): The temperature profiles of Fig. 14(b) have significantly higher 
uncertainty due to slight condensation in the jet. All other data have uncertainty mentioned in the text. 
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